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ABSTRACT 

We report the discovery by the HATSouth network of HATS-7b, a transiting Super-Neptune with a 
mass of 0.120zb0.012Mj, a radius of 0.563 ^q q 34 Rj, and an orbital period of 3.1853 days. The host star 
is a moderately bright {V = 13.340 zb 0.010 mag, Ks = 10.976 zb 0.026 mag) K dwarf star with a mass 
of 0.849zb0.027MQ, a radius of 0.815to;o35 -^0^ ^ metallicity of [Fe/H]= -b0.250zb0.080. The star 

is photometrically quiet to within the precision of the HATSonth measnrements and has low RV jitter. 
HATS-7b is the second smallest radius planet discovered by a wide-held ground-based transit survey, 
and one of only a handful of Neptune-size planets with mass and radius determined to 10% precision. 
Theoretical modeling of HATS-7b yields a hydrogen-helium fraction of 18 zb 4% (rock-iron core and 
H 2 -He envelope), or 9 zb 4% (ice core and H 2 -He envelope), i.e. it has a composition broadly similar 
to that of Uranus and Neptune, and very different from that of Saturn, which has 75% of its mass in 
H 2 -He. Based on a sample of transiting exoplanets with accurately (<20%) determined parameters, 
we establish approximate power-law relations for the envelopes of the mass-density distribution of 
exoplanets. HATS-7b, which, together with the recently discovered HATS-8b, is one of the hrst two 
transiting super-Neptnnes discovered in the Southern sky, is a prime target for additional follow-up 
observations with Southern hemisphere facilities to characterize the atmospheres of Super-Neptunes 
(which we define as objects with mass greater than that of Neptune, and smaller than halfway between 
that of Neptune and Saturn, i.e. 0.054Mj < Mp < 0.18Mj). 

Subject headings: planetary systems — stars: individnal (HATS-7) — techniques: spectroscopic, 
photometric 
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1. INTRODUCTION 

The HATSouth project (|Bakos et alJl201^ is the first 
global network of homogeneous, fully automated tele¬ 
scopes, capable of round-the-clock monitoring of a wide 
field-of-view (FOV) on the sky. 

The key scientific goal of HATSouth is to search for 
transiting extrasolar planets (TEPs), especially ones that 
are smaller and/or have longer periods than the hot 
Jupiters and Saturns, which have been the primary type 
of planet found by ground-based surveys to date. By cov- 
ering a wider total F OV than the NASA Kepler mission 
(jBorucki et al.l[2M^ . HATSouth monitors more bright 
stars per unit time than this mission, promising a greater 
yield of planets amenable to detailed characterization. 

The HATSouth telescopes are installed at Las Cam- 
panas Observatory (LCO) in Chile, at the H.E.S.S. site 
in Namibia, and at Siding Spring Observatory (SSO) in 
Australia. HATSouth was commissioned in late 2009, 
with regular operations beginning in 2011 after a ^1 year 
shakedown period, and has since collected 2.9 million 
science images at 4 minute cadence for some 9.5 million 
stars with r < 16 mag (2% per-point precision) in the 
Southern sky. As of May 1, 2015, the HATSouth tele¬ 
scopes have opened on 1788, 1528, and 1324 nights from 
LCO, HESS, and SSO, respectively. Based on weather 
statistics through 2015, the sites have averaged 8.24 hrs, 
7.54 hrs, and 5.24hrs of useful dark hours per 24 hr time 
period respectively. The longitudinal distribution of the 
three sites means that selected regions of the sky may be 
observed for long contiguous stretches of time, occasion¬ 
ally reaching 130 hours with interruptions shorter than 
30 minutes. 

Here we present the discovery of HATS-7b, which 
together with th e recently announced HATS-8b 
(|Bavliss et al.l 1201511 . is one of the first two transiting 
Super-Neptunes found by the HATSouth network, 
and one of only four such planets discovered by a 
ground-based wi de-held survey (t he other two be- 
ing HAT-P-llb. iBakos et al.l 1201(11 and HAT-P-26b, 
iHartman et al.ll2011bl l. By combining photometric ob¬ 
servations with high-resolution, high velocity-precision 
spectroscopy, we determine the mass and radius of 
HATS-7b to better than 10% accuracy. HATS-7b is one 
of only nine planets with M < 0.18 Mj for which both 
the mass and radius have been determined to this level 
of accurac}{3. With V = 13.340 ± 0.010 mag, HATS-7 is 
the 7th brightest star known to host a transiting planet 
with M < 0.18 Mj, and for which the mass has been 
accurately (with 10% acc uracy) determined . And of 
these planets, only G J 436 (jGillon et al.l l200 7tl. G.T 3470 
jBo^ls et al.l 1201211 and HAT-P-26 ( Hartman et al.l 
12011 HI haye deeper transits, making HATS-7b one of 
the best known low-mass transiting planets for detailed 
follow-up studies, e.g. transmission spectroscopy for 
probing its atmospheric composition. 

In the following section we describe the observations 
utilized to discover and characterize HATS-7b. In §[3]we 
discuss our analysis of the data to rule out blend scenar¬ 
ios and determine the system parameters. We discuss 
the results in § H) 

Based on the NASA exoplanet archive 
http://exoplaiietarchive.ipac.caltech.edu accessed 2015 
March 25. 
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Fig. 1.— Unbinned instrumental r band light curve of HATS- 
7 folded with the period P = 3.1853150 days resulting from the 
global fit described in Section [3] The solid line shows the best-fit 
transit model (see Section EJ. In the lower panel we zoom-in on 
the transit; the dark filled points here show the light curve binned 
in phase using a bin-size of 0.002. 


2. OBSERVATIONS 

2.1. Photometry 

2.1.1. Photometric detection 

The star HATS-7 (2MASS 13552567-2112276) was ob¬ 
served by the HATSouth wide-held telescope network 
between UT 2011 March 24 and UT 2011 August 19. 
These observations are summarized in Table [T] Observa¬ 
tions were made from LCO in Chile, the H.E.S.S. site in 
Namibia, and SSO in Australia. The data were reduced 
to trend-hltered light curves following iPenev et al.l ()2013[ 1 
and making use of the Trend Filtering Algorithm (TFA; 
iKovacs et al.l 1200^ . We searched the light curves ex¬ 
tracted from the images for periodic box-shaped transits 
using the Box Least Squares (BLS; IKovacs et al.l 1200^ 
algorithm, and identihed a P = 3.1853 day signal with 
a depth of 5.1mmag in the light curve of HATS-7 (Fig¬ 
ure [U the data are available in Table [5]) . 

We searched the HATSouth light curve of HATS-7 for 
other periodic transit signals which would indicate the 
presence of additional transiting planets in the system, 
but found no evidence for any such signals. We also 
searched for continuous periodic variability, due for ex¬ 
ample to the presence of star-spots on HATS-7. No sig¬ 
nificant periodic signals are present in the TFA-processed 
light curve, and we place an upper limit of 1 mmag on the 
peak-to-peak amplitude of any such variability. We also 
checked for continuous periodic variability in the light 
curve before applying the TFA filtering, and in this case 
find that the maximum signal present {P = 39 days) has 
a peak-to-peak amplitude of 2.7mmag, and is consistent 
with the level of low-frequency systematic noise present 
in the pre-TFA light curves of other bright stars in the 
field of HATS-7. 


2.1.2. Photometric follow-up 

Photometric follow-up observations of HATS-7 were 
performed using the 1-m telescopes in the Las Cumbres 
Observatory Global Telescope network (LCOGT) and 
the GROND instrument on the MPG 2.2 m telescope at 
La Silla Observatory (LSO) in Chile. These observations 
are summarized in Table [TJ Table [5] provides the light 
curve data, while the light curves are compared to our 
best-fit model in Figure [2j 


































HATS-7b 


3 


TABLE 1 

Summary of Photometric Observations of HATS-7. 


Facility/Field ^ 

Date Range 

Number of Points Median Cadence 
(seconds) 

Filter 

Precision ^ 
mmag 

HS-2/G568 

2011 Mar-2011 Aug 

5018 

290 

Sloan r 

7.1 

HS-4/G568 

2011 Jul-2011 Aug 

839 

301 

Sloan r 

6.9 

HS-6/G568 

2011 Mar-2011 May 

1627 

295 

Sloan r 

6.3 

LCOGTlm-l-SBIG 

2014 Jun 12 

110 

136 

Sloan i' 

2.0 

LCOGTlm-|-sinistro 

2014 Jun 16 

49 

167 

Sloan i' 

1.3 

GROND 

2014 Jul 20 

83 

140 

Sloan g 

1.2 

GROND 

2014 Jul 20 

83 

140 

Sloan r 

0.9 

GROND 

2014 Jul 20 

80 

140 

Sloan i 

1.0 

GROND 

2014 Jul 20 

80 

140 

Sloan z 

1.0 


^ For the HATSouth observations we list the HATSouth (HS) instrument used to perform the observations and 
the pointing on the sky. HS-2 is located at Las Campanas Observatory in Chile, HS-4 at the H.E.S.S. gamma- 
ray telescope site in Namibia, and HS-6 at Siding Spring Observatory in Australia. Field G568 is one of 838 
discrete pointings used to tile the sky for the HATNet and HATSouth projects. This particular field is centered 
at R.A. 14 hr and Dec. —22.5°. 

^ The r.m.s. scatter of the residuals from our best fit transit model for each light curve at the cadence indicated 
in the Table. 


Tw o transits were obse rved using the LCOGT 1 m net¬ 
work (jBrown et al.l 1201^ . The first, on UT 2014 June 
12, was observed using a 4Kx4K SBIG STX-16803 cam¬ 
era, with 0'/23 pixels, at the South African Astronomical 
Observatory (SAAO) station. The second transit, on UT 
2014 June 16, was observed using the Sinistro camera at 
Cerro Tololo Inter-American Observatory (CTIO) sta¬ 
tion in Chile, which utilizes a 4Kx4K Fairchild CCD-486 
back-side illuminated detector with 0'.'39 pixels. For both 
transits we used a Sloan i' filter. Because no apparent 
stellar neighbor is visible on archival images, or through 
our own observations, we defocused the telescope to re¬ 
duce systematic errors in the photometry. Standard 
CCD calibrations were performed and light curves were 
extracted using standard aperture photometry routines. 

A single transit was observed on UT 2014 July 20 us¬ 
ing the CROND multi-filter instrument, providing si¬ 
multaneous observations in the Sloan g, r, i, and z 
band-passes (|Creiner et al.l 1200811. The da t a wer e re- 
duced to light curves follow ing IPenev et al.l (1201311 and 
IMohler-Fischer et al.l (|2013ll . 

2.2. Spectroscopy 


Spectroscopic follow-up observations of HATS-7 were 
carrie d out with WiFeS o n the ANU 2.3 m telescope at 


SSO (iDooita et al.l 

200 

1), EEROS on the MFC 2.2m 

(IKaufer & Pasauinil 

1991 

j)' CORALIE on the Euler 1.2 m 

at LSO (lOueloz et al.l 2 

00 If), and HIRES on the Keck- 


110 m t elescope at Mann a Kea Observatory (MKO) in 
Hawaii (|Vogt et al.lll994ll . 


The WiFeS observations were carried out as part of 
our reconnaissance to rule ou t false positives, an d re- 
duced and a.nalyz ed following iBavliss et al.l (|2013ll and 
I Zhou et all (1201411 . We obtained three observations at 
a resolution of i? = A/AA = 7000 to check for radial 
velocity variations in excess of 5kms“^, and found that 
the observations were consistent with no variation at this 
level. A fourth spectrum at a resolution of i? = 3000 was 
obtained to provide an initial spectral classification and 
estimate of the surface gravity. Based on this observa¬ 
tion we found that HATS-7 is a K dwarf star with an 
estimated effective temperature of Teff* = 4700 ± 300 K 
and a surface gravity of log g* = 4.4 ± 0.3. 

We obtained three R = 48,000 resolution FEROS 



Time from transit center (days) 

Fig. 2. — Left: Unbinned follow-up transit light curves of HATS- 
7. The dates, filters and instruments used for each event are in¬ 
dicated. The light curves have been detrended using the EPD 
process. Curves after the first are shifted for clarity. Our best fit 
is shown by the solid lines. 

spectra and five R = 60,000 resolution CORALIE 
spectra, whi c h wer e reduced and analyzed following 
[Jordan et al.l (1201411 . The RVs measured from these 
spectra had a r.m.s. scatter of 20ms“^ and were con¬ 
sistent with no variation. Based on these observations 
we concluded that HATS-7b is a probable Neptune-mass 
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TABLE 2 

Differential photometry of HATS-7 


BJD 

(2 400 000+) 

Mag'^ 

o-Mag Mag(orig)‘’ 

Filter 

Instrument 

55784.52625 

-0.00947 

0.00553 

r 

HS 

55685.78275 

0.00387 

0.00490 

r 

HS 

55682.59768 

0.00210 

0.00476 

r 

HS 

55666.67146 

0.00587 

0.00552 

r 

HS 

55762.23228 

-0.00517 

0.00421 

r 

HS 

55784.52962 

0.01251 

0.00600 

r 

HS 

55650.74655 

0.01291 

0.00467 

r 

HS 

55749.49145 

0.00899 

0.00508 

r 

HS 

55685.78731 

0.00285 

0.00494 

r 

HS 

55682.60232 

-0.00122 

0.00467 

r 

HS 


Note. — This table is available in a machine-readable form in the online jour¬ 
nal. A portion is shown here for guidance regarding its form and content. The 
data are also available on the HATSouth website at http://www.hatsouth.org 
^ The out-of-transit level has been subtracted. For the HATSouth light curve 
(rows with “HS” in the Instrument column), these magnitudes have been de¬ 
trended using the EPD and TEA procedures prior to fitting a transit model 
to the light curve. Primarily as a result of this detrending, but also due to 
blending from neighbors, the apparent HATSouth transit depth is somewhat 
shallower than that of the true depth in the Sloan r filter (the apparent depth is 
85% that of the true depth). For the follow-up light curves (rows with an Instru¬ 
ment other than “HS”) these magnitudes have been detrended with the EPD 
procedure, carried out simultaneously with the transit fit (the transit shape is 
preserved in this process). 

^ Raw magnitude values without application of the EPD procedure. This is 
only reported for the follow-up light curves. 


planet, and continued observing it with the larger aper¬ 
ture Keck-I telescope. The stellar atmospheric param¬ 
eters estimated from the FEROS and CORALIE spec¬ 
tra are consistent with those measured from WiFeS. We 
also determined that the star is slowly rotating (with 
usinz < 5kms“^), and, based on inspecting the cross¬ 
correlation functions (CCFs), found no evidence for ad¬ 
ditional stellar components in the spectra. 

We observed HATS-7 with the HIRES spectrometer 
(Vogt et al. 1994) on the 10-m Keck-1 telescope using 
standard practices o f the California Planet Survey (CPS; 
iHoward et al.ll20lO) . Our 10 observations through a gas 
cell of molecular iodine (I 2 ) lasted 25 minutes each us¬ 
ing the “C2” decker (14 x 0.86 arcsec slit). In addi¬ 
tion, we observed HATS-7 without the I 2 cell using the 
“B3” decker (14 x 0.57 arcsec slit) to record a tem¬ 
plate spectrum for the RV analysis and for measuring 
high-precision stellar atmos pheric parameters. We com¬ 
puted relative RVs using the lButler et al.l (1199611 method 
to model the I 2 x star spectrum. RV errors were es¬ 
timated from the uncertainty on the mean of 700 spec¬ 
tral segments (each spanning ^2 A) that were separately 
analyzed for each observation. We also measured spec- 
tral line bisector s pans (BSs) from these data following 
[Torres et all (j2007[ ). The final RVs and BSs are provided 
in Tableland are displayed in Figure O The RVs vary in 
phase with the photometric ephemeris and with a semi¬ 
amplitude of AT = 18.4 ± 1.9ms“^, while the BSs are 
consistent with no variation and have an r.m.s. scatter 
of 8.4ms“^. 

Based on the Keck/HIRES spectra, HATS-7 is a quiet 
K dwarf, with barely detectable Calcium HK activity 
{{S) = 0.27, logi?(jK = -4.800; Tab. H]). Also, HATS- 
7 has no significant RV jitter (Tab. [S]). It is simi- 
lar to other very low activity K dwarfs, like HAT-P-26 
(|Hartman et alfcoildl . 



Fig. 3. — Top panel: High-precision RV measurements for 
HATS-7 from Keck/HIRES, together with our best-fit circular or¬ 
bit model. Zero phase corresponds to the time of mid-transit. The 
center-of-mass velocity has been subtracted. Second panel: Veloc¬ 
ity O—C residuals from the best-fit model. The error bars for each 
instrument include the jitter which is determined in the fit. Third 
panel: Bisector spans (BS), with the mean value subtracted. Note 
the different vertical scales of the panels. 


3. ANALYSIS 

We analyzed the photometric and spectroscopic ob- 
servations of HATS-7 to determine the parameters of 
the system using the standard procedures develop ed for 
HATNet and HATSouth (see I Bakos et al.l l201(il with 
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TABLE 3 

Relative radial velocities and bisector span measurements of 

HATS-7. 


BJD 

(2 456 000-I-) 

RV''- 

(ms-l) 

(ms“^) 

BS 

(ms-l) 

O-BS 

Phase 

Instrument 

826.78970 

21.15 

2.07 

4.9 

2.5 

0.709 

Keck 

826.80787 

17.17 

2.21 

-3.3 

6.0 

0.715 

Keck 

827.78004 

-3.77 

2.46 

-3.2 

3.3 

0.020 

Keck 

827.79805 

-2.93 

2.47 

-8.6 

5.5 

0.026 

Keck 

827.81761'' 



13.2 

10.2 

0.032 

Keck 

828.78151 

-16.51 

2.22 

-4.2 

4.4 

0.334 

Keck 

828.79940 

-14.90 

2.35 

5.5 

7.5 

0.340 

Keck 

829.77933 

5.93 

2.67 

7.7 

9.7 

0.648 

Keck 

829.79741 

19.86 

2.62 

-12.2 

9.4 

0.653 

Keck 

830.77141 

3.32 

2.19 

-8.6 

6.1 

0.959 

Keck 

830.78899 

6.17 

1.88 

8.9 

6.9 

0.965 

Keck 


^ The zero-point of these velocities is arbitrary. An overall offset 7rei is 
fitted to the Keck/HIRES velocities (Section O, and has been subtracted. 

^ Internal errors excluding the component of astrophysical/instrumental 
jitter considered in Section fSl 

^ This HIRES observation was taken without the iodine cell to be used as 
a template. The RV is not measured for this observations, but the BS value 
is measured. 


0.5 



4.0 


Effective temperature [K] 

Fig. 4.— Comparison between the measured values of Teff* and 
p* (from SPC applied to the HIRES R-free template spectrum, and 
from our modeling of the light curves and RV data, respectively), 
and the model isochrones from [Yi et al.l 11200l|j . The best-fit 
values (dark filled circle), and approximate Icr and 2a confidence 
ellipsoids are shown. The values from our initial SPC iteration are 
shown with the open triangle. The isochrones are shown for 
ages of 0.2 Gyr, and 1.0 to 14.0 Gyr in 1 Gyr increments. 

modifications described bv iHartman et al.]l2012H . 

High-precision stellar atmospheric parameters were 
measured from the Keck/HIRES template spectrum us- 
ing the Spectral Param eter Classification program (SPC; 
iBuchhave et al.l[2Q12all . The resulting Teg* and [Fe/H] 
measurements were combined with the stellar density p* 
determined through our joint light curve and RV curve 
analysis, to determine the stellar mass, radius, age, lu¬ 
minosity, and other physic al parameters, by comparison 
with the Yonsei-Yale (Y^; lYi et al.l l2001f) stellar evolu¬ 
tion models (see Figure 0]). This provided a revised es¬ 
timate of log p* which was fixed in a second iteration of 
SPC. Our final adopted stellar parameters are listed in 
Table 01 We find that the star HATS-7 has a mass of 
0.849 ± 0.027Mq, a radius of 0.815tg;Q35 i?©, and is at 
a reddening-corrected distance of 261^12 pc- Note that 
the error on the distance does not take into account 
the systematic errors due to uncertainties in the stellar 
isochrones. HATS-7 is a main sequence star, and has a 
poorly constrained age due to the limited stellar evolu¬ 


tion expected within the age of the universe at this stellar 
mass. 

We also carried out a joint analysis of the HIRES 
RVs (fit using a Keplerian orbit) and the HATSouth, 
LCOGT 1 m, and G ROND light curves (fit using a 
IMandel fc Agoll [20021 transit model wit h fixed quadr atic 
limb darkening coefficients taken from IClaretl l2004f) to 
measure the stellar density, as well as the orbital and 
planetary parameters. This analysis makes use of a dif¬ 
ferential evol ution Markov Ch ain Monte Carlo procedure 
(DEMCMC; Iter Braa^l2006fl to estimate the posterior 
parameter distributions, which we use to determine the 
median parameter values and their Icr uncertainties. We 
als o varied the jitte r as a f ree parameter in the fit follow¬ 
ing iHartman et al.l (|2014f) , where we use the empirical 
distribution of the Keck/HIRES RV jit ter values for l ow- 
activity cool main sequence stars from iWrighd (|2005ll to 
place a prior on the jitter. The results are listed in Ta¬ 
ble 0] We find that the planet HATS-7b has a mass of 
0.120±0.012Mj, andaradiusof0.563t)j;)!gi?j. We also 
find that the observations are consistent with a circular 
orbit. When the eccentricity is allowed to vary in the fit, 
we find e = 0.065 ± 0.053 and a 95% confidence upper- 
limit of e < 0.170. The parameters listed in Table[S]were 
determined assuming a fixed circular orbit. 

In order to rule out the possibility that HATS-7 is 
a blended stellar eclipsing binary system, we carried 
out a blend analysis of the photometric data following 
IHartman et al.l (|2012f) . We find that a model consisting 
of single star with a transiting planet fits the data better 
than any of the blended stellar eclipsing binary models 
that we tested, but there are some blend models which 
cannot be rejected with greater than 5cr confidence based 
on the photometry alone. In these models the primary 
star in the background eclipsing binary has a distance 
modulus that is no more than 2.1 mag greater than that 
of a brighter foreground star. We simulated the cross¬ 
correlation functions, RVs and BSs that would be mea¬ 
sured with HIRES/Keck-I for these blend scenarios and 
found that in almost all cases the object would be easily 
identified as a composite stellar system (the CCFs would 
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be double-peaked, or the BSs and/or RVs would vary by 
greater than 100ms“^). We searched the spectrum for 
th e spectral si g nature of a secondary star as described 
in iKolbl et al.l (|2015ll . and found no companions down 
to 1% of the brightness of the primary at separations in 
RV greater than ±10kms“^. There is a small region of 
parameter space where the BS and RV scatter for the 
blended system would be comparable to the measured 
scatters. These are blends where the background binary 
has a distance modulus that is between 1.95 mag and 
2.1 mag greater than that of the foreground star, and 
where the RV of the foreground star differs by at least 
50kms“^ from the systemic RV of the binary. Consid¬ 
ering, however, that these blends can still be rejected 
with 4(7 confidence based on the photometry, and that 
the simulated RVs do not show the sinusoidal variation 
in orbital phase that is observed, we conclude that these 
blend scenarios do not provide a good description of the 
observations, and that HATS-7 is a transiting planet sys¬ 
tem. 

While we can rule out the possibility that HATS-7 is 
a blended stellar eclipsing binary system, we cannot rule 
out the possibility that HATS-7 is a transiting planet 
system with a fainter unresolved stellar companion. We 
note that this is often the case for many known tran¬ 
siting exoplanet systems. We simulated this scenario as 
well, and find that companions with masses up to that of 
the planet host cannot be ruled out. The highest spatial 
resolution images we have available are from the HIRES 
guide camera based on which we can rule out compan¬ 
ions with AV < 4 at a separation greater than 1.5". 
High mass companions would also need to have been ob¬ 
served near conjunction for the system not to have been 
detected as a spectroscopic binary. Higher spatial reso¬ 
lution imaging and/or additional RV monitoring would 
be needed to put tighter constraints on any stellar com¬ 
panions. If HATS-7 does have a stellar companion, the 
radius and mass of the planet HATS-7b would be some¬ 
what larger than what we infer here. 


4. DISCUSSION 

With its mass of Mp = 0.120 ± 0.012 Mj, HATS-7b 
is one of the very few transiting and well characterized 
super-Neptunes known to date. We adopt a nomencla¬ 
ture, whereby planets with 0.18 Mj < Mp < 0.3 Mj, 
i.e. those more massive than halfway between Neptune 
and Saturn, are referred to as sub-Saturns (esp. those 
toward the lower end of this mass range), whereas plan¬ 
ets with 0.054 Mj < Mp < 0.18 Mj as super-Neptunes 
(esp. those toward the higher mass end). This nomen¬ 
clature is more or less c onsistent with previou s lit¬ 
erature, e.g. HAT-P-18b (jHartman et al.M2nilal l is a 
called a sub-Saturn ( Mr, = 0.197Mj), while Kepler- 


101b (|Bonomo et al.l 1201411 is called a super-Neptune 
(Mp = 0.16 Mj). The mass vs. radius of low mass 
(Mp < 0.2 Mj) transiting planets with accurately de¬ 
termined masses (here we use <20% uncertainty) is 
shown in Fig. [S] It is noteworthy that only a hand¬ 
ful of super-Neptunes are known (with accurate param¬ 
eters), and the regime between 0.1 Mj and 0.18 Mj is 
occupied by only 5 s uch objects. The se are Kepler- 
9c (K=12. 34 mag; [Torres et al.l 1201111 . Kepler-lOlb 
(K=12.0; iBonorno et al.l 201411. the rece ntly discov¬ 
ered HATS-8b fK=12.6: iBavliss et al.ll2015ll , Kepler-35b 


TABLE 4 

Stellar Parameters for HATS-7 


Parameter 

Value 

Source 

Identifying Information 



R.A. (h:m:s) 

13^55^25.688 

2MASS 

Dec. (d:m:s) 

-21°12'27.7" 

2MASS 

R.A.p.m. (mas/yr) 

-43.3 ± 1.4 

2MASS 

Dec.p.m. (mas/yr) 

-37.2 ± 1.9 

2MASS 

GSC ID 

GSG 6148-00422 

GSC 

2MASS ID 

13552567-2112276 

2MASS 

Spectroscopic properties 



Teff* (K) . 

4985 ± 50 

SPC ^ 

Spectral type. 

K2 

SPC 

[Fe/H]. 

+0.250 ± 0.080 

SPC 

ilsinf (kms ^) ... 

0.50 + 0.50 

SPC 

7rv (kms-l) . 

9.435 ± 0.017 

EEROS 

(s> . 

0.27 

Keck/HIRES 

log ^hk . 

-4.800 

Keck/HIRES 

Photometric properties 



B (mag). 

14.359 ± 0.020 

APASS 

V (mag). 

13.340 + 0.010 

APASS 

9 (mag). 

13.825 + 0.010 

APASS 

r (mag). 

13.013 + 0.030 

APASS 

i (mag). 

12.690 + 0.030 

APASS 

J (mag). 

11.528 + 0.024 

2MASS 

H (mag). 

11.085 + 0.024 

2MASS 

Ks (mag). 

10.976 + 0.026 

2MASS 

Derived properties 



MMMq) . 

0.849 + 0.027 

Y2+p*+SPC 

K,(Rq) . 

n R1 e:+U.u4y 

U.OlO_Q Q35 

Y2+p*+SPC 

P* (cgs) . 

2.22 + 0.33 

y2+p*+spc 

logs* (cgs). 

4.545 + 0.049 

Y2+p*+SPC 

LMLq) . 

0.368^^38 

Y2+p*+SPC 

My (mag). 

6.10 + 0.14 

Y2+p*+SPC 

Mx (magjESO) 

3.94 + 0.11 

Y2+p*+SPC 

Age (Gyr). 

7.8 + 5.0 

Y2+p*+SPC 

Ay (mag) <=. 

0.186 + 0.058 

Y2+p*+SPC 

Distance (pc). 

257/j® 

Y2+p*+SPC 


^ SPC = “Stellar Parameter Classification” method based on 
cross-correlating high-resolution spectra against synthetic templates 
iBuchhave et ^12012^ . These parameters rely primarily on SPC, 
but have a small dependence also on the iterative analysis incor¬ 
porating the isochrone search and global modeling of the data, as 
described in the text. 

^ Isochrones-|-p*-hSPC = Based on the isochrones et al.l 
nsni), the stellar density used as a luminosity indicator, and the 
SPC results. 

^ Total V band extinction to the star determined by comparing the 
catalog broad-band photometry listed in the table to the expected 
magnitudes from the Isochrones-|-p*-|-SPC model for the star. We 
use the lCardelli et al.l II1989D extinction law. 


(K=13.9; I Welsh et akllMll . and HATS-7b (K=10.97), 
the subject of this paper. Among these planet host¬ 
ing stars, HATS-7 is by far the brightest, which will 
greatly facilitate detailed studies. All low-mass transit¬ 
ing planets shown in Fig. [5] were found by either Kepler 
and Corot (space-born surveys), HATNet and HATSouth 
(ground-based wide-field surveys), th e pointed ground- 
based survey MEarth (G J 1214b; lirwin et al.l l2009t 
iCharbonneau et abl I2009H . or they were found by the 
photometric follow-up of planetary systems that were 
first discovered by RVs. 

Well characterized planets and host stars in this regime 
are important in exploring the correlations between var¬ 
ious system parameters. One example is the host star 
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Fig. 5.— Mass—radius diagram of super-Neptunes {Mp < O.lSMj) and super-Earths with accurately measured masses and radii (< 20% 
uncertainties). Color indicates equilibrium temperature (with the palette of R,G,B = 2000, 1500, 1000 K), while size is inversely proportional 
to surface gravity. HATS-7b is also marked with a large box. Isodensity lines for p = 0.2, 0.5, 1.27 (Uranus), 1.64 (Neptune), 5.5 (Earth) 
gcm“^ are shown with dashed lines. Neptune is marked with a blue triangle. Abbreviations are: K: Kepler, H: HAT, HS: HATSouth, C: 
Corot. 


metallicity - planet occurrence rate correlation, which 
has been known for a long time (jFischer fc Valentill2005[l 
to exist for Jupiter size planets, and which correlation 
wa s found to be weakened for small, < 47?® planets 
by jBuchhave_et_^n ( 2012bll . However, recent analysis 
by I Wang fc Fischer! (|2015ll shows that the planet occur¬ 
rence rate-metallicity correlation is universal, including 
for terrestrial planets. Using an inc reased sample and 
new analysis, iBuchhave et al.l (l2014fl find a correlation 
between host star metallicity and planet radius, in the 
sense that the average metallicity of the host star in¬ 
creases with planet size. These authors also find three 
populations of exoplanets (rocky, gas-dwarf, giant), as 
based on the metallicity of the stars. The discovery of 
HATS-7b around a metal rich ([Fe/H]= -1-0.250 ± 0.080) 
star, and the apparent lack of a hot Jupiter in this system 
will help in refining such future analyses. 

The mass regime of super Neptunes is important 
in studying the transition from ice giants to gas 
giants, where, in t he core accretion scenario (e.g. 
iMordasini et al.l[2015[l . rapid accumulation of a gaseous 
envelope is expected to start. This transition is well 
demonstrated in Fig. [51 which plots the bulk density 
of transiting exoplanets as a function of their masses 
(again, only for those with reliable parameters deter¬ 
mined to better than 20% precision). Mean density first 
decreases with increasing planetary mass, as the plan¬ 
etary mass increases from super-Earths ('^ 0.01 Mj) to 
Neptunes (^ 0.1 Mj). Then, in the super-Neptune/sub- 
Saturn regime, the trend reverses, as matter becomes 


increasingly degenerate under high pressure, and density 
starts to increase. 

With the increasing number of well characterized tran¬ 
siting exoplanets, we are gradually mapping out the pa¬ 
rameter space these planets occupy. While it is custom¬ 
ary to explore the relations between various parameters 
(e.g. mass, radius, equilibrium temperature), here we ex¬ 
amine the boundaries spanned in the planetary mass — 
density domain. The approximate envelopes of the dis¬ 
tribution are: 


Mp < 0.02 Mj 


Pp >14(Mp/Mj)0-i™ 
Pp <25(Mp/Mj)0-i78 


i.e. for rocky planets without significant amount of 
volatiles. 


Mp < 0.4 Mj 


Pp > 3/4 • 10-4 (Mp/Mj)-2 
Pp < 3/5 (Mp/Mj)-3/4 


i.e. for Neptunes and Saturns, and 


Mp > 0.4 Mj 


Pp > 1/5 (Mp/Mj)3/2 
Pp <5/2(Mp/Mj)i-i 


for ice and gas giants {pp measured in gcm”^ for all). 
Here the pp oc (Mti/ Mj)°- 4^^ relation f or the smallest 
planets is taken from iSotin et al.l (|2007[ 1. and is based 
on models. The other relations are approximate, and are 
not physically motivated, though it is possible that phys¬ 
ical explanations for them do exist. Between earths and 
super-earths, the mean density only increases slightly. 
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TABLE 5 

Parameters for the transiting planet HATS-7b. 


Parameter Value ^ 


Light curve parameters 

P (days) . 

Tc (BJD) b . 

Ti 4 (days) ^ . 

7 i 2 = T 34 (days) .. 
a/Ri, . 

C/R*‘= . 

Hp/i?* . 

. 

b = a cos i / Ri, . 

i (deg) . 


3.1853150 ±0.0000054 
2456528.29697 ± 0.00058 
0.0958 ±0.0012 
0.00735 ± 0.00099 

lU.O»_Q 
22.62 ±0.16 
0.0711 ± 0.0019 

U.i4o_Q Qg 5 
+ 0.11 


0.39_g 43 
87.92 ±0.75 


Limb-darkening coefficients 


ci,i (linear term) . 0.4140 

C 2 ,* (quadratic term) .. 0.2464 

ci,r . 0.5519 

C 2 ,r . 0.1984 


RV parameters 

K (ms-l) . 

e ® . 

RV jitter (ms“l) ^ . 

Planetary parameters 

Mp (Mj) . 

Pp {Pj ) . 

C{Mp,Rp) s . 

Pp (gcm-3) . 

logflp (cgs) . 

a (AU) . 

Teg (K) h . 

0 ' . 

(F) (10^ergs“l cm“^) ' 


18.4 ±1.9 
< 0.170 
3.0 ± 1.1 


0.120 ±0.012 


0.563 


-1-0.046 

-0.034 


- 0.02 


0.83 ±0.18 
2.968 ± 0.076 
0.04012 ± 0.00043 


1084 ± 32 


0.0200 ± 0.0024 


3.12 


+0.45 

-0.32 


^ The adopted parameters assume a circular orbit. Based 
on the Bayesian evidence ratio we find that this model is 
strongly preferred over a model in which the eccentricity is 
allowed to vary in the fit. For each parameter we give the 
median value and 68.3% (Icr) confidence intervals from the 
posterior distribution. 

° Reported times are in Barycentric Julian Date calculated 
directly from UTC, without correction for leap seconds. Tc'. 
Reference epoch of mid transit that minimizes the corre¬ 
lation with the orbital period. T 14 : total transit duration, 
time between first to last contact; T 12 = T 34 : ingress/egress 
time, time between first and second, or third and fourth con¬ 
tact. 

^ Reciprocal of the half duration of the transit used as a 
jump parameter in our MCMC analysis in place of a/R*. It 
is related to a/R* by the exp ression C/^* — Q'/P*(27r(l + 
e sin u;))/(P-\/l — b'^y/l — e^) JBakos et al.|[2010D . 

^ Values for a quadratic law, adopted from the tabulations 
byEI are 3 t2nnli according to the spectroscopic (SPC) pa¬ 
rameters listed in Table [T] 

® The 95% confidence upper-limit on the eccentricity from 
a model in which the eccentricity is allowed to vary in the 
fit. 

^ Error term, either astrophysical or instrumental in origin, 
added in quadrature to the formal RV errors for the listed 
instrument. This term is varied in the fit assuming a prior 
inversely proportional to the jitter. 

® Correlation coefficient between the planetary mass 
Mp and radius Rp determined from the parameter pos¬ 
terior distribution via C{Mp, Rp) — {{Mp — {Mp)){Rp — 
{Rp}))/{(7Mp<yRp)), where (•) is the expectation value op¬ 
erator, and (Tx is the standard deviation of parameter x. 

^ Planet equilibrium temperature averaged over the orbit, 
calculated assuming a Bond albedo of zero, and that fiux is 
re-radiated from the full planet surface. 

^ The Safronov number is given by © — ■^(V esc/V)rb)^ = 
(a/Rp)(Mp/M*) fsee lHansen fc Barman||2007D . 

^ Incoming fiux per unit surface area, averaged over the 
orbit. 


Once Mp > 0.02 Mj ( 6 M 0 ) is reached, the bulk density 
decreases as a result of retaining more volatiles, and the 
ability to form a more extended atmosphere. Finally, 
at around Mp > 0.2 Mj 60 M 0 ), the compression of 
matter gets significant to result in an increasing mean 
density. 

The lower limit for the density of gaseous planets 
(dashed line in Fig.jBJ oc is related to the inflation 

mechanism of planets, which has not been fully solved. 
It is noteworthy that this relation holds over 2 orders of 
magnitude in mass, from Mp w 0.5 Mj to 50Mj. The 
upper limit for the density of rocky planets (dotted line, 

oc Mp is likely to suffer from the bias against discov¬ 
ering compact rocky/iron planets, due to their smaller 
transit signatures. Also, while there are close to 250 well 
characterized transiting exoplanets, considering the di¬ 
mensionality of the parameter space (e.g. Mp, Rp, Mi,, 
P, [Fe/H], Teff*, age), this is still a very small number 
to understand their distribution. Consequently, signifi¬ 
cant biases in the observed distribution exist. Probably 
the most important such bias is that the majority of the 
planets plotted in Fig. [0] have short periods. The confir¬ 
mation and characterization of long period planets takes 
much more time, and these systems are primarily dis¬ 
covered by Kepler, so the host stars are typically faint, 
further hindering follow-up studies. 

The bulk composition of HATS-7b may be estimated 
by fitting theoretical structural and evolutionary mod¬ 
els to the observed mass and radius of 38.0 ± 3.9 M^ 
( 0.120 ± 0 . 012 Mj) and i ?0 (0.563i;j;^^^ i?j). 

We compute a range of models for planets of HATS-7b’s 
mass f or various possible param eters using the methodol¬ 
ogy of lHowe fc Burro-^ (|2015f) . If a two-layer partition¬ 
ing with a rock-iron core and a hydrogen-helium enve¬ 
lope is assumed, the best fit models suggest a core mass 
of 31±4 M 0 and an envelope mass of 7±1.5 M 0 , that 
is, a hydrogen-helium fraction of 18±4%. This result is 
consistent within uncertainty across a range of metallic- 
ity and whether or not evaporative mass loss is incor¬ 
porated in the model. More quantitatively, increasing 
the metallicity from lx solar to 10 x solar increases the 
radius at 7.8 Gyr by about 0.2 i ?0 (0.018i?j), approxi¬ 
mately the same amount as decreasing the core mass by 
1 M 0 . This is also of the same order as the uncertainty 
introduced by the uncertainty in age. Models with sig¬ 
nificantly smaller cores, <25 M 0 , are ruled out because 
they under no circumstances shrink to the observed ra¬ 
dius of the planet over the age of the system. For com¬ 
parison, final radii for the 10 x solar models with core 
masses of 0, 5, 10, 15, and 20 M 0 are 0.95, 0.87, 0.79, 
0.71, and 0.65 Rj, respectively (all larger than the ob¬ 
served radius of 0.56 i?j), while a bare core of 38 M^ has 
a radius of 0.202 i?j if composed of rock-iron and 0.296 
Rj if composed of ice. In the case of a pure ice core, the 
best fit is a core mass of 34.5±4 M 0 and an envelope 
mass of 3.5±1.5 M^, that is, a hydrogen-helium fraction 
of 9± 4%, similar to the composition of Uranus and Nep¬ 
tune l|Helled et al.ll20HT) . Fig.[3shows radius versus time 
for representative evolutionary models with core masses 
of 28-34 M 0 , encompassing the full range of core masses 
consistent with observations. 

While the mean density of exoplanets is an important 
constraint on their bulk composition, detailed studies 
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Fig. 6. — Planetary mean density as a function of planetary mass for planets with mass measured to better than 20% precision. HATS-7b 
is also marked with a large box. The size of the points scales with planetary radius, while the color indicates equilibrium temperature (with 
the palette of R,G,B = 2000, 1500, 1000 K). It is clearly visible that the less dense Jupiter-mass planets are predominantly ‘hot’ (red). The 
four thin lines indicate approximate bo undaries of the c urrent transiting extrasolar planet population in this space. The significant outlier 
on the bottom left s ide is Kepler-87c llOfir et al.ll201'Jl . and the other significant outlier at Mp ^ 1 Mj and p ^ 3gcm~^ is WASP-59b 
IIHebrard et al.|[2013l') . Solar system planets are marked with blue triangles. 



Fig. 7.— Radius versus time for representative evolutionary 
models plotted against the observed radius and age estimates 
for HATS-7b. Models are labeled with their corresponding core 
masses. We include models with core masses from 30 to 34 
with a H 2 -He atmosphere with solar metallicity, from 28 to 34 
with an atmosphere with lOx solar metallicity, and four best-fit 
models with 30x solar metallicity. In all cases, we fit the models 
to a (final) total mass of 38 . 

of the planetary atmospheres can provide a direct mea¬ 
surement of the composition of the planet. One impor¬ 
tant way to detect such compositions is via transmis¬ 
sion spectroscopy. In Fig. |5]we show super-Neptunes to 
super-Earths, including the sub-Saturn HAT-P-18b, with 
the expected transmission signatures plotted against the 
planetary masses. To derive an approximate measure 


of the transmission signatures, we approximated the H 
scale height as H = kBTeq/{ngp), where fcs is the Boltz¬ 
mann constant, T^q is the equilibrium temperature of 
the planets (assuming black bodies, and full redistribu¬ 
tion of heat), fi is the mean molecular weight, and gp 
is the planetary surface gravity. To make the compar¬ 
ison simpler (and because of lacking further informa¬ 
tion on their atmospheric composition), we used p, = 2 
(i.e. H 2 ) for all planets, even though the mean molecu¬ 
lar weight of the atmospheres will be larger, and thus 
our estimates are somewhat optimistic. We then ap¬ 
proximated the fractional contribution of transmission 
through the planet’s atm osphere to the tot al stellar flux 
as (5 = 5 X 2RpH/R^ (iPerrvmanI I2014D . and multi¬ 
plied this with the total estimated K-band flux of the 
star. Atmospheric transmission analysis has been car¬ 
ried out for five of these planets, typically those in the 
top left side of Fig. |8] (those with low mass and high 
expected transmission signal). These are, in order of 
incre asing planetary mass: GJ-1214b (iKreidberg et ^ 
[Ml. HD 97658b (ll^itson et all l2014bD. G.T 3470b 


(lEhrenreich et al.l2011. GJ 436b (iKnntson et al.ll2014l 


and HAT-P-llb ( Fraine et al.l 12011 . With the excep¬ 

tion of HAT-P-llb, all showed essentially featureless 
transmission spectra, indicating hazes, clouds, or at¬ 
mospheres with high molecular weight. The only Nep¬ 
tune mass planet with features in its transmission spec¬ 
trum is HAT-P-llb, where the signature of water was 
detected ([Frame et al.l l20ll . Among the super Nep- 
tunes, HATS-7b is the most promising target for de¬ 
tecting its atmosphere. Also, such measurements will be 
facilitated by the quiet host star; high resolution spec¬ 
troscopy of HATS-7 shows very low RV jitter, low v sini, 
and no chromospheric activity. Altogether, it is scien- 
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tifically compelling to pursue transmission spectroscopy 
for the recently discovered super Neptunes HATS-7b and 
HATS-8b, and the sub-Saturn HAT-P-18b, where rela¬ 
tively large signals are expected. Perhaps a trend from 
featureless spectra to feature-rich transmission (observed 
for certain hot Jupiters) will be observed with increasing 
planetary mass through the super-Neptune regime. 
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Fig. 8. — Approximate detectability of a planet’s atmosphere in transmission as a function of planetary mass for super-Neptunes 
{Mp < O.lSMj) and super-Earths . In calculating the scale height of the planet atmosphere we assume the molecular weight of pure 
molecular hydrogen, and the equilibrium temperature and surface gravity of the planet determined as from our analysis (Tab. [5ll. The 
fractional contribution of the transmission signal to the star light is approximated as 5 = 5 x 2RpH/R^ ll Perrvmanll 201 411 . We then calculate 
the K-band flux of the star, and multiply with 5, to come up with an approximate measure of the transmission signal. We stress that this 
quantity does not take into account the detailed (expected) spectrum of the planetary atmosphere, but is an order of magnitude estimate 
of the signature. The size of the filled circles scales with the radius of the planet (arbitrary scale), and the radius of the open circles scales 
with the scale height of the atmosphere. Small plus symbols denote planets with uncertain mass or radius measurements (error >20%). 
Abbreviations are: K: Kepler, H: HAT, HS: HATSouth, C: Corot. 
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